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Abstract. In the focus of our attention is the asymptotic properties of the 
sequence of convex hulls which arise as a result of a peeling procedure applied 
to the convex hull generated by a Poisson point process. Processes of the 
considered type are tightly connected with empirical point processes and stable 
random vectors. Results are given about the limit shape of the convex hulls in 
the case of a discrete spectral measure. We give some numerical experiments 
to illustrate the peeling procedure for a larger class of Poisson point processes. 
Keywords : Control measure ; convex hull ; limiting shape ; peeling ; Poisson 
point processes ; stable vectors. 



1. Introduction 

Consider a Poisson point process (p. p. p.) tt = 'i^a,v with points scattered over 
W^. Identify \ {0} with x S'^'^ where S'^'^ is the unit sphere. Assume that 
the intensity measure of this process /i is of the form 

(1) 12 = 9 XV 

where 

(1) 6 is the absolutely continuous measure on M.]_ determined by the density 
function 

^W-ar-"-i, r>0, 
dA 

A is the Lebesgue measure in M^, a > is a parameter while 

(2) ly, called the spectral density, is a bounded measure on the cr— algebra Bgd-i 
of the Borel subsets of S'^~^. Without loss of generality we assume that 
i'{S'^^^) = 1. We denote by Si, the support of v. 

The representation ([T]) means that for any Borel A C and E C S'^^^ 



|a; |a;| e A, £ E^ = e{A)iy{E). 



where ex = |a;|^"'"a;, for all x £ M"*-. 

We assume that the Poisson point process tTq^^ is non-unilateral. It means that 
ly is supported by a set Si, C S"^^^ such that the cone cone{S^) generated by Si, 
coincides with M.'^. 
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Let 5(0, r) denote the ball of a radius r centred at the origm. It is easily seen 
that for any S > 

/i(S(0,5)) = oo, fi{R'^\B{0,S)) < oo. 

It implies that with probability 1 in any neighbourhood of the origin there are 
infinitely many points of tt while 7r(E'* \ B{0,6)) is finite. 

The interest to the point processes controlled by ([T]) is explained by the following 
facts. 

Let ^(1), ^(1), ^(") be independent copies of a random vector ^ such that 
the function P{|^| > r} regularly varies as r — > oo with the exponent —a and the 
measures {i'r)r defined by 

i^riE) = P{e^ eE : \^\>r}, Ee S^.-i, 

weakly converge to v on B'^{0, r) for any r > 0. 

Consider the empirical point process /3„ generated by ^^"^ or, 

more precisely, by the random set 

,l"U{6-ie(^),...,&-ie(«)} = UW,...,|(")} 

where 

6„ = inf {r : nP{e^£E\ \^\ > r) < l} . 

It is easily seen tlicit the point process [S^i weakly converges to tt — tTq, (see e.g. 
[5], Prop. 3.21). Thus, each weak limit of a sequence of the empirical 

processes. 

It can be easily established that na^i, admits the following representation 

oo 

(2) ^"^- = E'^{r-^^°e,} 

fe=l 

where 

k 

9 Tk = ^^7fej the sequence {'jk)keN is a sequence of i.i.d. random variable 

1=1 

with common exponential distribution with mean equal to one. 

• (efe)feGN is a sequence of i.i.d. with common distribution v 

• {lk)keN and {ek)k(EN are supposed to be independent 

(see e.g. |2j). It is worth recalling that the point processes tTq.i^ naturally arise 
within the framework of the theory of stable distributions. For example, if a G 
(0,1) then the series ( = J^xUie-rr ^^''^ converges a.s.. Furthermore, Q has the 
d— dimensional stable distribution ( see e.g. [5], [B], [5]). 

In the focus of our interest is the sequence of convex hulls that arise from the 
peeling procedure introduced in [T]. In what follows by C(A), yl C ffi.'^, we denote 
the convex hull generated by A. Let C be a convex set. By ext C we denote the 
set of the extreme points of C. If C is a convex polyhedron then ext C is the finite 
set of its vertices. 

It convenient to start with the binomial process (<;^"''). Let cj""* = C(<j^"''). If 
the measure P^, the distribution of ^, has no atoms then a.s. C^"'' is a polyhedron 
and, furthermore, 

n a cl"^ = ext 
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Define 

4") = ^(") \ ext = 

then 

and so on. Obviously, tlie sequence of tlie so-built non empty convex liulls C^"' , C2"'' , • • • , C^"'' , 
is finite and its length is random. 

Definition 1. We say that the underlying distribution and corresponding to 
it spectral measure v are non-unilateral ij the minimal closed cone containing S^, 
coincides with . 

If the underlying distribution is non-unilateral then for any fixed k 

e int C^"' , inf I a; e 9(7^"' } > a.s. 

for n sufficiently large. In [IJ, it was shown that if < a < 2 and v is non- 
unilateral then (C^"\ Cf\ C^"'), as n — >■ 00, converge in distribution to 
(Co, Ci, Cfe) for any fixed k. Consequenlty, the sequence #{extC^"''} is 
bounded in probability as n tends to infinity . 

in) 

In order to learn how Cj. relates to Ck when k — kn ^ 00 we need, first, to learn 

(n) 

how C'k behaves as fc — cx). It should be noted that can be regarded as the 
multi-dimensional analogue of the order statistics. So, the asymptotic properties of 



cl"-* are of great interest from the view-point of mathematical statistics. 



We generalize now the construction of the peeling sequence to infinite set. 
Let <r = <ri denote the set of points of tt or, in other words, let <j support the random 
measure tt. We may apply to ^ the same peeling procedure as in the case of the 
finite set As a result we obtain the sequence of sets <ji, ^2, • ■ • , the sequence of 
their convex hulls Ci, C2, ■ ■ ■ and the sets of extreme points ext Ck, k = 1,2, . . . . 
Furthermore, a.s. <;k+i = <?fc \ ext Ck- If is non-unilateral then € Ck a.s. for all 
k. Furthermore, Ck is a.s. a polyhedron and 

<?fc n a Cfc = ext Ck. 

Intuitively, we expect that the asymptotic behaviour of C„ is rather regular. It 
is convenient to state our basic conjecture in the following way: 
Denote 

(3) Cn = Pn^Cn, where p„ = max \x\. 

If I' is non-unilateral then there exists a non-random set 
C such that 

lim dniCnjC) — a.s. 

Definition 2. C (if it exists) is called the limit shape of the sequence Cn. 

It is easy to show that if such C exists then it is certainly non-random. Indeed, 
labelling the points of <^ in the descending order of their distances from the origin 
we obtain a sequence a;'-^' , x*-^' , x'^^^ , . . . such that a.s. 

> > > ... . 
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It is worth noting that the joint distribution of jx'-"'^-'!, |a;^^^|, . . . , is absolutely 

continuous with the density of the form 

(4) p„(ri, ra, . . . ,r„) = a"(rir2 • • • r-„)"""^e-'^(^''^')''r°I{^^>^^>...>^,j. 
Let (ry, e), (ryi, e'^-*), {r]2, e'^-*),... be i.i.d. with common distribution 

P{77 > r, eeE}^ e"V(£;). 
According to ([2]), we have 

(5) {x^'nT=i ^ + • ■ ■ + 

which implies that the event {lim„_>oo Cn exists} belongs to the cr-algebra I of 
the events invariant with respect to all finite permutations of the random vectors 
(7?i, e(i)), {ri2, By the Hewitt-Savage zero-one law I is trivial. Since 

the limit set C — lim„^oo C'n is X-measurable we conclude that C is constant with 
probability 1. 

Now we give a first example where the existence of the limit shape is proved. 

Example 1. Let Si, consist ofd+1 unit vectors e'^^\ . . . ^e'^'^^'^^ such that cone{e'-^\ . . . ,6^'*"''"'^-') 
coincides with M.'^. iVa.v is decomposed on d+ \ one- dimensional independent p.p.p. 

oftheform ^x^'-* — l^^.*'' je^'^^ . Since v is non-unilateral the points x'"f!\ i — 1, 2, . . . , d-l- 

1, serve as vertices ofC'k, k = 1,2,... . Moreover, \x^^\ = (i'i)^^"('7i + • • ■ + ?/fe)"'^^" 
and pnn^^" — > i+ = maxKi^d+i a.s.. Then the limit shape C is the convex 

polyhedron with vertices u^*-* = {ti/t'^)e^^^ and ti = {viY'/'^ , i = 1, 2, . . . , d + 1. 

If #(>5'y) > d -\- I then the situation becomes much more complicated. Theorem 
[T] and [2] proved below deal with a case where a non-unilateral v is supported by a 
finite number of unit vectors. 

Intuitively, we expect that, say, in case of ly uniformly distributed over S'^~^ 
the unit ball arises as the limit shape. However, it is not easy at all to prove this 
formally. The authors tried to verify the credibility of this conjecture using the 
Monte Carlo simulation. Obviously, the representation ([2]) provides a basis for such 
a simulation. The results of simulation presented below make this conjecture very 
credible. 

It should be emphasised that the basic goal of the present paper is to draw 
attention to new and interesting problems of stochastic geometry. So far, little or 
nothing is known about the peels no matter what point process they concern. 

The paper is organised as follows. In Section 2, we obtained a partial result 
on the limit shape of the convex hulls Ck{Tra,u) when the spectral measure of the 
process tTq is atomic. Section 4 contains some numerical experiments. 

2. Almost sure convergence of the peeling 

In this section we assume that the spectral measure v of the process tTq, is 
atomic, i.e. it is supported by a finite number of the points e'^-*, . . . , e''-* belonging 
to the unit sphere S'^~^. 

Furthermore, it is also assumed that cone{e'-^\ . . . , e'-'^} = R'^. Denote by j/^ = 
^({e^*-'}), i = 1,2,. .. ,1, the atoms of v. 
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It implies that the considered point process is a superposition of the one-dimensional 
independent Poisson point processes defined on the rays 

^{x\ ie^*', i > 0}, 

If a Borel set A C Ci then 

(6) fi{A) = v,a [ r-"-Mr. 

J A 

Definition 3. Let A = {a^^^ . . . , a^™)} be a finite set where a^*^ G M'' for i = 
1, . . . , m, and m > d + 1. The set A is extreme if 

ext C{A) = A. 

Theorem 1. Let Ck{T^a.v) be the k-th convex hull of the Poisson point process TTa^v 

Denote by Coo the convex hull generated by A = {h'l^"e^^\ . . . , vj^"' e^''"'}. 
If A is extreme then as fc oo 



(7) dH (fc^/"Cfe(^„^,), Coo) = O ( 1 a.s. 



where Unfc — Inlnfc. The polyhedron Coo determines the limit shape of the convex 

hulls Ck{T^a,v). 

Remark 1. If a is uniformly distributed over its support in the sense that Vi = l~^ , 
then the total number of the vertices of Ceo equals I. Furthermore, they lie on the 
sphere of the radius /^^/". Loosely speaking, the convex hulls Ck{iTa,v) are getting 
round as k oo. 

If the condition A is extreme is omitted, we can state the following result : 

Theorem 2. Let Ck{T^a,u) be the k-th convex hull of TTa.u- Denote by Coo the 
convex hull generated by A — {v\^'^e^^\ . . . , h'f^°'e^''^}. Then as k ^ oo 

Coo I —5* a.s. 



(8) dH (A:i/"Cfe(7r„,,), Coo) 



Remark 2. Let f be a continuous homogeneous functional of a degree 7 defined 
on convex sets. From Theorem\^ we get 

fiCkij^a^v)) /(Coo) a.s. 

In particular, if f{A) is the surface Lebesgue measure, i.e. f{A) = X'^^^{dA), 
then 

(9) /(Cfc(^„,,))- 



k^ 



But if f{A) = X^iA), then 

f{Co 



(10) /(Cfc(^„,,)) 



6 



YU. DAVYDOVS A.V. NAGAEV ^ AND A. PHILIPPE^ 



3. Proof 

3.1. Auxiliary lemmas. Let ?7i,r/2, ... be i.i.d. random variables with the stan- 
dard exponential distribution, so that a = E771 = Var?7i — 1. Define the sums 

(11) r„ = 771 H 

By the law of the iterated logarithm there exists an a.s. finite random variable k 
with values in N such that for n > k 



(12) |n-ir„ - 1| < 2W— a.s.. 

V n 

Consider a function h{z) = z^^/". If |z — 1| < 1/2 then 
\h{z) ~ h{l)\ <La\z-l\, Lo,< 00. 

Let n' = min |n| 2^J^ < 1/2| . If n > niax(n', k) then 



\h{n-'r^)-h{l)\<2L^J^^ 



and, therefore, for n > k = k(w) 



(13) 



V Unn 



We call the configuration any countable set of points from M*^ such that for any 
i5 > there are a finite number of points belonging to the set that lie outside 
the ball {a; | |a;| < 5}. So the point is the limit point of any configuration. 
We call a configuration non-unilateral if all the convex hulls, Ck = Ck{(^), k = 
1,2,..., generated by <; contain as an interior point. It is evident that under the 
conditions of Theorem [T] the random measure ttq.^^ is almost surely supported by a 
non-unilateral configuration 

Denote by int(<r) the set of the interior points of <j, i.e. 

int(<;) = {a; I X e x 9Ci(<^)}. 
Lemma 1. Let cri, <J2 G be such that q C then for all k Cz N, 

(14) Cfe(^i) c Cfc(^2). 

Proof. It is trivial that Ci(<;i) C Ci(<;2)- 

Note that if x is an interior point of Ci(iri), i.e. x G int(Ci(<ji)), then x is also 
an interior point of Ci(<;2), therefore 

int(Ci(^i))Cint(Ci(<j2)) 

and this implies that 

C2(<?i) = Ci(int(Ci(^i))) C Ci(int(Ci(^2))) = ^2(^2) 

By induction, the lemma is proved . □ 

Lemma 2. LetK, he the set of non-unilateral configurations such that no d+l points 
lie on the same hyperplane. Let ^, ^' € /C he such that ^' C <j and ^{<;'\(,) — m < 00. 
Then we have, for all k £ N, 

(15) Cfc+„,(OcCfe(<j')cCfc(?) 
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Proof. Since e /C and is the only limit point of both configurations all 

Ck{<i), k = 1,2, . . . , are polyhedrons. Note that for all k,l > 1 

(16) Cfe+i(0 = Ci(int(Cfe(OnO) 
and 

(17) Cfc+,(0 = Cfc(int(Ci(^)nO). 

First, let m = 1. Note that the inclusion Ci(<;') C Ci(<j) follows directly from 
the relation C Denote {a} = \ Consider two possible cases a ^ Ci(<;') and 
a e Ci{';') one after another. 

Let a ^ Ci{(^'). In this case a £ (9Ci(<j) i.e. Ci(<j') ^ Ci('?)- It implies that 
int(<j) C Utilising (|16p under k — 1 yields C2(^') C Ci(<;). Since the inclusion 
Ci(?') C Ci(?) is obvious we conclude that ([T5|) holds for fc = m = 1. 

Further, let us make use of the induction by k. Assume that ([T5|) holds for m — 1 
and all fc < n and show that then it holds for to = 1 and k = n + 1. By the induction 
assumption we have 

(18) C„+i(^) C C„(^') c C„(^). 
Since 

int(C„(<j')nO =int(c„(<j')n^) 
we obtain, taking into account (|16p . 

c„+i(0 = c^i(int(a.(OnO)- 

From the right hand side inclusion of ((T8)) . it follows that C„+i(^') C C„+i(^). 
Further, from the left hand side inclusion of (fT8|) we conclude that 

int(C„+i(0 n = int(C„(<;') n ^'). 

Applying ([171) yields C„+2(<r) C C„+i(^')- Thus, ^ holds for fc n + 1 and to = 1, 
i.e. the case a ^ Ci(<;') is exhausted. 

If a G Ci (<;'), then there exists an integer uq such that 

C„(?')-C„(0, n = l,2,...,no 

Furthermore, a ^ C„q+i (<;'). Obviously, the relations (|15p are trivial for to = 1 and 
n = l,2,...,no. Hence, it remains to apply the above argument to the configu- 
rations C„(,+i(<;) n <; and C„o+i(?') n Thus, the lemma is proved for all k and 
TO = 1. 

Now, let m > 1, i.e. \ = {oi, . . . , Um}- 
Consider the configurations 

'^0 = = \ {Ol}, 

"^2 = ? \ {01,02}, ■ • ■ , 

W = ^ \ {oi, . . . ,a„i} = 

Note that the neighbouring configurations differ by a single point. So, one may 
apply pS)) . Applying it yields 

Cfc+m('?) C Cfc+m_i(<;i) C Cfc+m_2(T2) C . . . C Cki'^rn) — Ck{<^') C Cfc(<;). 

The lemma is proved. □ 
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3.2. Proof of Theorem [H 

Lemma 3. Let e^*^ for i — 1, . . . ,1 with I > d + 1, be unit vectors such that 

cone {e(i\...,eW}==M''. 

If A ~ {i'l^'^e^^\ . . . , vj^^e^'^} is extreme then there exists r > and e depending 
only on A and on the dimension d such that the set {rii>l^" e^^\ . . . yrivj^" e^''^} is 
extreme for all (ri, . . . , r„) such that \ri/r — 1| < £, i — 1, . . . ^l. 

Proof. The proof of this lemma is evident. □ 

Proof of Theorem [7J Let us label the points lying on the ray Ci in the descend- 
ing order of their norms. So, we have the sequence x^^\x2\--- such that a.s. 
|a;^*^| > \x2^\ > ■ ■ ■ . Obviously, the sequences {xn\n £ N}, 1 < i < I, are jointly 
independent. Furthermore, from ([6]) 

where r„ is defined as in ([TT|) . 

Let e > 0. According to there exists no = rio{uj) such that for all z = 1, . . . , Z 
and all n > hq 



(19) 
and 



Let the configuration is formed by the points Xn\ n > no, i ~ 1, . . . , / i.e. 

c' - I I-Tt^*) t''^ \ 

Consider for all A; < 1 

A+ = {(no + k- l)-i/"(l + e,)i^!/"e(i), . . . , (no + fc - l)-i/"(l + efc)^.y"e«} 
and 

A- = {(no + k- l)-i/"(l - efe)i^!/"e(i), . . . , (no + fc - l)-i/"(l - Sk)i^f^"e('^} 
where 



'lln(fc + no - 1) 



fc + no — 1 

By virtue of (HH) the points . . . , xill hit the layer C{A'^) \ C{A];). Then by 
Lemma 13] the convex hull Ci(<;') is the polyhedron and 

extCi(0={xW,...,a;W}. 

Similarly, the set 

{^io + l' ■ • ■ '^ii + li ^i^t) \ C(^2 ) 

and, therefore, it is extreme, i.e. 

ext C{{x[^^_^_l, . . . , = {xl^g^^i, ■ ■ ■ , ^nl+i}- 

It is evident that 

extC2(0 = {4?+i,---,4?+i}- 
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Continuing in this way we obtain at the k-th convex hull Cki'i') such that 

ext Cu{,') = {41.-1' • • -^^no+u-i) C C{At) \ C(A-). 
The last inclusion implies that 

dH ((fc + no - Coo) < £fc, 

where, we remind, Coo is the convex hull generated hy A — {v\^°'e^^\ . . . , v^^°'e^''^}. 
From (Uni it follows that 

Ck+,nW) C Ck+7n{T^a,i^) C Cfe (<;'), with m ^ {uq - 1)1. 

Therefore, 

ext Cfc+™(7r„,^) C CiA+) \ C{A^^^,), with m' = (no - !)(/ - 1). 
So, for all sufhciently large k 

dH ((fc + m)i/"Cfe+„(7r„,,), Coo) < 2sk. 
Since m is fixed the theorem follows. □ 

3.3. Proof of Theorem [2j Let e be an arbitrary positive real. Hereafter, we 
denote ^^'^ the e-neighbourhood of a set A, 

A^''> = {x : d{x,A) < e}. 



Let Ai = An dC{A) = {^-j^^e^J), j e J}, the set Ai is extreme. From the 
process TTa.i^ , we construct a new p. p. p. tti obtained by deleting all the points on 
the rays Cj = {x\ x = te^'^\ t > 0}, j £ J. By Lemma [1] we have for all n € N 

(20) C„(^i) C C„(^o,.) 

Moreover, Ai is extreme and Coo — C{A) ~ C{Ai), thus Theorem [1] ensures the 
convergence 

(21) dff (ni/"C„ (tti), Coo) -^Oa.s. 
From (PO)) and (PT|) . it exists ni € N such that for all n > ni 

(22) Coo C ni/"C„(7ri)(^) C ni/"C„(7r„,,)('). 
It is easy to see that there exists i £ I such that the set. 

je J; ;z>y"e«, z G {1, ...,/} \ J} 
is extreme and satisfies the following relation 

(23) CooCC(A2)cC(:). 

From -Ka.u we construct a second p. p. p. 7r2 by adding the independent point 
processes {iiiji^j verifying the following conditions 

• the {TTi)i^j are independent of iTa^u '■, 

• for each i e J the spectral measure of Hi is supported by Ci and the intensity 
measure is fj,i{A) = (i>i — i^i)a /^r^"^^d r. 
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According to Theorem [T] we have 

(24) dff(ni/"C„(7r2),C7(A2)^0a.s. 
and using Lemma [1] we get, for ah n € N, 

(25) C„(7r„,.) C Cn{^2)- 

From and (UHl), it exists 712 G N such that for ah n > n2 

(26) ni/"C„(^„,,) C ni/"C„(^2) C C(A2)(^) C C^^'\ 
According to ([221) and (|26p . for aU n > max(7ii, 712), we have 

n^^^CniiTa,.) CCg'^ and Coo C 7ii/"C„(^„^,)(2'). 
By definition of dn , this means 

),Coo) < 2e, 

and we get ([5]). 

4. Simulation and conjectures 

We investigate using some simulations the hmit shape and the asymptotic be- 
haviour of basic functionals in the case of continuous spectral measure. Hereafter 
we consider the example of the uniform distribution as spectral measure. 

The point processes {x*-^^ , j S N} are simulated using the representation ([5]). 
Let Ci,„ be the convex hull generated by the first n points x^^\ x'^^\ . . . jX*^"^ and 
Kn.i = min \x\. Since the points of the simulated p. p. p. are ordered by their 

distances from the origin, it is evident that 

Ci^n' = Ci with n' = min{n : k„^i > 

This fact is used to construct the successive convex hulls {Ck)k&h 

Figures [T] gives an impression about the behaviour of the peels. The observed 
closeness of the peels to the unit circle also support our conjecture about the exis- 
tence of the limit shape that is expected to be a circle. 

It is of great interest to get impression about a possible behaviour of such basic 
functionals of the convex polygons Ck, fc = 1, 2, . . . , as the perimeter C, the area A 
and the total number of vertices M. It seems evident that C{Ck) and A(Ck) tend 
to zero as A; — >■ 00. Intuitively, we expect that N{Ck) — )■ 00 as fc — > 00. Figure 
[2Jtop] represents the logarithm of those functionals as function of log(fc), calculated 
on simulated p. p. p. for different values of a. The observed closeness of the points 
to straight lines makes it reasonable to expect that in a sense 

(27) C(Ck) - k--", A{Ck) - k-^\ N(Ck) x fc^" 

with 7i and 7^ depending on a whereas it seems that 7„ does not dependent on a. 

The next step consists in estimating the exponents and possibly the dependence 
on a. Using independent replications of p. p. p., we estimate the three exponents 
defined in ([27]) for different values of a. Figure [2] [bottom] represents the logarithm 
of the estimated exponents versus log(a). For the three cases, the linear approx- 
imation seems reasonable. According to the estimated coefficients of the straight 
lines (see the equations in the caption of Figure [2]) , it looks very credible that the 
true values are 

3 3 1 

(28) 7; = — , 7^ = - and 7„ = -. 

Za a 2. 
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After the k-th iterative step of the peehiig procedure, the number of deleted 
points should be the order of X]^=i-^(^j) ^'^'^ 

Pk — max Ixl X I fc2 

Moreover we can expect that dn {p^^Ck{T^a.u), C'oo) converges to zero. 

Using the arguments of Remark [21 this convergence would lead to C{Ck) ^ Pk^ 
^-3/(2a) A{Ck) ^ Pk ^ fc~^/". These convergence rates are in agreement with 
the estimated values of ai a,d aa obtained in (|28p . 



25 50 




T "h 1 ] 1 |J T 1 T 1 
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Figure 1. [dotted line] the normalized simulated shapes of C25, 
C'so; C'loo and C150 in the case a — 3/2 and the spectral measure 
is uniform, [solid line[ the unit circle. 
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Figure 2. [TOP] Log- log representation of the values of 
£{Ck), A{Ck) and Af{Ck) as function of k. The functionals are cal- 
culated on simulated p. p. p. for different values of a = 0.5; 1.0; 1.5 
and the uniform distribution as spectral measure. [BOTTOM] The 
estimated values of the logarithm of exponents defined in ((27)) ver- 
sus ln(a) and the best linear fittings In 7a — —0.97 In a -I- In 2.95, 
In 7, = -0.971na + lnl.48 and ln7„ = 0.06 Ina + lnO.48. Three 
exponents are estimated for each a on 1000 independent replica- 
tions. 
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